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We introduce a natural and simple to implement regular- 
ization scheme of the Hartree-Fock-Bogoliubov (HFB) equa- 
tions with zero range pairing interaction. This renormaliza- 
tion scheme proves to be equivalent to a simple energy cut-off 
with a position dependent running coupling constant. 

PACS numbers: 21.60. Jz, 21.60.-n, 21.30.Fe 



More than forty years after the pioneering work of 
Bohr, Mottelson and Pines |^ there is no need to re- 
iterate again the relevance of pairing correlations in nu- 
clei. It is well established that nuclei are s-wave "super- 
conductors" in the so called weak coupling limit, when 
the pairing gap is much smaller then the fermi energy 
A ^ Ei? = h'^kp/2m {kp is the fermi wave vector). In 
this limit one can show that the rms radius of the Cooper 
pair (in infinite matter) significantly exceeds the interpar- 
ticle separation, h'^kp/mA ^ 1/kp and the radius of 
the nucleon-nucleon interaction as well. As in the case 
of the deuteron, the details of the two-particle interac- 
tion at distances smaller or comparable with the inter- 
action radius should be irrelevant and the bulk Cooper 
pair properties should be described basically by a single 
constant, derivable from a suitably chosen zero range in- 
teraction model. (We shall not address here corrections 
beyond the leading order, such as effective range effects.) 
One encounters typically no insurmountable difficulties 
in introducing a local HF (or Kohn-Sham) Hamiltonian 
h{r) ^j^- If one can adopt the approximation of a zero 
range two-body interaction in the pairing channel as well 
then the HFB equations become 

[h{r) - ti]u,{r) + A{r)v,{r) = E.u^ir), (1) 
A*(r)M,(r) - [h{r) - ^i]v,{r) = E,v,{r). (2) 

Here Ui{r) and Vi{r) are the quasi-particle wave func- 
tions, A(r) = is the local pairing field, fi is the 
chemical potential, Egs is the ground state energy of the 
system and v{r) is the anomalous density. In all the 
formulas presented here we shall not display the spin de- 
grees of freedom. If one takes at face value Eqs. (|l|,||) one 
can show that the diagonal part of the anomalous den- 
sity matrix v(r, r) diverges, since when |ri — ^0 the 
anomalous density v{ri,r2) has the singular behavior 



(3) 



and the local pairing field A{R) cannot be defined [|-0. 

In metals this type of singularity does not play a no- 
ticeable role, because the summation over the single- 
particle states is cut-off at energies of the order of the 
Debye energy lod £f- The single-particle density 
of states is essentially constant in an energy window of 
width 0{ijJd) cf and the expression for the anomalous 
density has only an infrared logarithmic divergence. This 
logarithmic divergence is due to states near the Fermi 
surface and has nothing to do with the ultraviolet diver- 
gence due to states faraway from the Fermi surface, which 
leads to the 1/|ti — singularity discussed here. The 
infrared divergence leads to the notorious non-analytical 
dependence of the gap on the coupling constant, namely 
A = wri exp(— l/yTV), where V is the strength of the 
interaction and N is the single-particle density of states 
at the Fermi energy ep. 

In nuclei and especially in very dilute fermionic atomic 
systems, where kpr^ <ti 1 and tq is the radius of the inter- 
action, there is effectively no well defined cut-off and one 
needs to regularize the theory. A finite range interaction 
will provide a natural cut-off at single-particle energies 
of the order of Sc It? /mr'^, when the fast spatial oscil- 
lations of the quasi-particle wave functions Ui{r),Vi{r) 
will render the nonlocal pairing field A(ri,r2) ineffec- 
tive. Even though the presence of a finite range of the 
interaction in the pairing channel formally removes the 
ultraviolet divergence of the gap, it is very difficult to 
come to terms with the fact that a cut-off at an energy 
of the order of Ti' /mr^ could be the responsible for the 
definition of the gap both in the case of regular nuclei 
and very dilute nuclear matter as well. The character- 
istic depth of the nucleon-nucleon interaction potential, 
which is of the order of Yi' lmr\, being the largest energy 
in the system, can be effectively considered to be infinite 
in the case of dilute systems. A well defined theoretical 
scheme for the calculation of a local pairing field, should 
lead to a converged result when only single-particle states 
near the Fermi surface are taken into account. 

Most of the calculational schemes suggested so far for 
infinite systems reduce to replacing a zero range poten- 
tial by a low energy expansion of the vacuum two-body 
scattering amplitude The traditional approach 

in the calculations of finite nuclei consists however in in- 
troducing a simple energy cut-off, while the pairing field 
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is computed by the means of a pseudo-zero-range in- 
teraction. In this approach the effective range of the 
interaction is obviously determined by the value of the 
energy cut-off and the two-body coupling constant in 
the pairing channel is chosen accordingly Such a 

pure phenomenological approach lacks a solid theoretical 
underpinning and always leaves the reader with a feel- 
ing that "the dirt has been swept under the rug". An- 
other solution favored by other practitioners is to use a 
finite range two-body interaction from the outset, such as 
Gogny interaction . Besides the fact that the ensuing 
HFB equations are much more difhcult to solve numeri- 
cally, such an approach also lacks the elegance and trans- 
parency of a local treatment and this seemingly simple 
recipe is indeed as phenomenological in spirit as the treat- 
ment based on a pseudo-zero-range interaction, with an 
explicit energy cut-off. Moreover, in spite of the feeble 
arguments often put forward in favor of a finite range 
interaction in HFB calculations, the only real argument 
is the fact that the pairing field would otherwise diverge, 
and there is no mean-field observable which would be no- 
ticeable different in the case of a finite range interaction. 

The only attempt to implement a consistent regular- 
ization scheme for finite systems that we are aware of is 
that of Ref . . In agreement with the analysis of Ref. H] 
the authors of Ref. conclude that in the case of a zero 
range two-body interaction the anomalous density has a 
l/|ri — singularity. The regularization schemes for 
infinite homogeneous systems amounts to subtracting a 
term proportional to l/k"^ in the gap equation in momen- 
tum representation which in coordinate representa- 
tion corresponds naturally to a l/|ri — term as well. 
Since the divergence in the anomalous density v{ri,r2) 
is due to large momenta and thus short distances, it is 
not surprising that the character of the divergence is not 
affected by the size of the system. Bruun et al. advo- 
cate the use of the following calculational procedure for 
the anomalous density. First of all one represents the 
anomalous density as [18| 



A(r 



v*{ri)ui{r2) + 



Ei>0 

Goirl,r2,^i), 



A(r) V;*(ri)V>,(r2) 
2 fJ.~ Si 



\h{r) - £i]ipi{r) = 0, 

[/i - /i(ri)]Go(ri,r2,/i) = 5{ri - r2), 



(4) 

(5) 
(6) 



where r = (ri 4- t'2)/2. One can easily justify this sub- 
traction scheme in infinite homogeneous matter, since 
<(ri)u,(r2) = ^r^{rM^{r2)/2^{e, - + A^. In 
the limit ri r2 the sum over single-particle states 
in Eq. (jj) is converging now and one has only to extract 
the regulated part of the propagator Go(ri, r2, /i), using 
the pseudo-potential approach H 



v*{r)u^{r) + 



A(r)VA*(r)?A.(r) 



Ei>0 



m 



Gq lim Go(ri,r2,Ai) -I- 21 \ 

ri-*r2 2nh ri — r2 



(7) 

(8) 



obtaining for the local pairing field 



A(r) 



E 



A(r)V'*(r)V'^(r) 



47r|a|?i^ A(r) 



-Bi>0 



2in~ Si) 



(9) 



where a is the two-particle scattering length (a < 0). 
The renormalization procedure and the extraction of the 
regulated part from various diverging quantities is com- 
pletely analogous to the familiar procedures in Quantum 
Field Theory, with the only difference that in this case 
everything is performed in coordinate space. One liter- 
ally "throws away" the diverging terms and retains the 
nonvanishing finite contributions. 

The approach suggested in Ref. Q has however two, 
related, problems and as is formulated is applicable for 
systems in a harmonic trap only and does not apply to 
atomic nuclei or other self-sustaining systems. First of 
all, after the divergence has been eliminated, the regu- 
lated expressions for the anomalous density and for the 
pairing gap, Eqs. (0,^, are defined entirely in terms of 
states in a certain neighborhood of the fermi level, since 
the corresponding sums converge rather quickly. Only 
when one can establish a one-to-one correspondence be- 
tween the HFB terms v*{r)ui{r) and the corresponding 
HF expressions AirJ^ip* {r)^pi{r) /2{fi — Si) it is clear how 
to evaluate Eqs. There is no one-to-one correspon- 

dence for self-sustaining systems |H,|l9| , where sufficiently 
deep bound hole states lie in the continuum and where 
often there is no one-to-one correspondence between the 
HF and HFB spectra around the fermi level. In the case 
of nuclei very close to the nucleon drip lines the HFB 
spectra are continuous essentially everywhere, while the 
HF spectra are not, and the one-to-one correspondence 
between HFB and HF is absent. The second and the 
most difficult aspect of the approach suggested in Ref. 
[Q however is the fact that it requires the determination 
of the regular part of the single-particle Green function 
G^'^^ (R, j-i) , for which there is so far no clear computa- 
tional scheme in the case of an arbitrary self-consistent 
field. These two problems are to a large extent related, 
as only the whole expressions (0j9|) are uniquely defined, 
but not each separate part. 

Our suggestion amounts to a simple to implement ap- 
proach. First of all we introduce an explicit energy cut- 
off Ec in evaluating the anomalous density. In this way 
we can evaluate separately the HFB and HF sums in 
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Eqs. (y,0j8|j9|) irrespective of the existence of the one-to- 
one correspondence discussed above. The final result is 
independent of Ec, if this is chosen appropriately. Sec- 
ondly, we remark that there is no compelling reason to 
use the exact HF single-particle wave functions, energies 
and propagator in Eqs. and in order to construct 

the regulator one can use a Thomas-Fermi approxima- 
tion for the relevant quantities. Since the divergence has 
an ultraviolet character, the Thomas-Fermi approxima- 
tion is particularly well suited . Thus we arrive at the 
following relations 

m exp{ikp{r)\ri — r2|) 



Go(ri,r2,Ai - U{r)) 



27rr 



2T:h \ri — r2 
,(r) := iy,{r) + 



ikp{r)m 
iA{r)kF{r)m 



0(|ri 



^2 



r2 



A(r) 
47r2 



fce(r) 



k^dk- 



1 



2 1,2 



h^k 

2m 



U{r) -\- 17 



A{r)mkc{r) 



2n^h' 



1 - 



Vc{r) 
h{r)-- 



Ei<Ec 



2KXr) 
v*{r)u,{r), 



kpir) . kc{r) + kpir) 
In 



kcir) - kpir) 



2m 



h'klir) 



2m 



h'k'pjr) 
2m 



Uir) 
hC/(r) - 
U{r), 



(10) 

(11) 
(12) 

(13) 

(14) 
(15) 
(16) 



where the cut-off energy Ec is chosen sufficiently far away 
from the Fermi level to insure that the rhs of Eqs. ( pT|jT^ ) 
has converged. As usual one has to take the limit 7 0+ 
at the end of the calculations. The local wave vector 
kpir) is real only in the physically allowed region of the 
fermi level, where the regularized part of the propaga- 
tor is imaginary. This imaginary part of the regularized 
propagator is, naturally, exactly canceled by the corre- 
sponding imaginary part of the momentum truncated 
propagator in Eq. (pT|). If the Fermi momentum be- 
comes imaginary (outside nuclei for example) one can 
easily show that Vregif) is still real. The pairing field 
has thus the simple expression 



A(r) 
1 

9eff{r) ~ 
mkc{r) 

~ 2TT^n^ 



'9eff{r)Vc{l') = -gVreg{r) 

1 
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1 



kp{r) fcc(r) + kp{r) 



(17) 



(18) 



2kc{r) kc{r)-kp{r)' 

where g = ATrTi^a/m. Surprisingly, these relations look 
very much like a simple position or density dependent 



renormalization of the coupling constant. For a typical 
nuclear potential which monotonically increases with the 
radial coordinate {dU{r)/dr > 0) one can easily show 
that dgeff{r)/dr > 0, thus the effective pairing interac- 
tion is stronger inside than outside nuclear matter (re- 
member g < 0). This is stark contrast with the behav- 
ior one would get using the popular energy cut-off of a 
g6{ri — r2) interaction, namely the vacuum renormal- 
ization scheme In this case the effective coupling 
constant is gvaci'i') = g /{^ — gmkc{r)/2TT^h^] and one can 
then easily show that dgyac{r) / dr < if dU{r)/dr < 0. 

It is instructive to apply this recipe to the case of infi- 
nite homogeneous matter. After a few simple manipula- 
tions one can show that the equation for the gap reads 



1 

kp 



dk 



2fcF|a| 



V(fc2 - fc^)2 + fc4 

, , 2fcc|fl| 



kF\a\ kc 
m — 



kp 



(19) 



where fcp = 2m A /h^. Using the methods described in 
Refs. P,p|-p^ one would not get the term with the log- 
function. The technical reason is that we used A/(si — /z) 
instead of A/e^ in Eqs. (|0J|,|ll|) respectively |fl|, which 
enhances the convergence of the corresponding sums or 
integrals discussed above. Parametrically we are allowed 
to make such a substitution as long as |fci^-a| ^ 1, other- 
wise one should consider effective range corrections and 
higher partial waves. Even though the momentum cut- 
off kc appears explicitly here, once this momentum cut- 
off is sufficiently large, there is no dependence of the gap 
on the cut-off momentum. 

When evaluating the total energy of the system one 
has to be careful and calculate the expression |l2[ 



E„ 



2m 



Tc{r) - A{r)uc{r) 



E^, 



pot J 



(20) 



where Epot is the usual HF potential energy con- 
tribution, since the kinetic energy density Tc(r) — 
2 S_B<_B I^^-e('')I^ diverges in a similar fashion as I'dr) 
Ec, but Egs does not. 

We have implemented this renormalization scheme 
for the pairing field for both selfconsistent and non- 
selfconsistent calculations of spherical nuclei. The nor- 
mal and anomalous densities were computed following 
the complex energy integration technique extensively 
used by Fayans and his collaborators |]l^. In order to 
illustrate the convergence properties we present in Fig. 1 
the neutron pairing field A(r) obtained as a solution of 
the Eqs. (|^,|^,|l^,^ for a range of cut-off energies Ec- 
The calculations were performed for a simple Woods- 
Saxon potential with fixed parameters corresponding to 
a ^^^Sn nucleus and for a fixed value of the chemi- 
cal potential /i = —0.1 MeV (essentially at the neutron 
drip line). The value of the bare coupling constant is 
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g = —200 MeV- fm^. The total energy converges equally 
fast with Ec- The reasons why convergence is achieved 
for Ec = C(e_F) and how one can improve on this aspect 
are discussed in Ref. 1231. 




FIG. 1. The neutron pairing field ( [l7| ) as a function of 
the radial coordinate and of the cut-off energy Ec- Upward 
various curves correspond to Ec = 20, 30, 35, 40, 45 and 
50 MeV respectively. On the scale of the figure the last two 
curves are indistinguishable. 

In conclusion, wc have presented a renormalization 
procedure for the HFB equations in the case of zero range 
pairing interaction, which is easy to implement for any 
type of finite or infinite systems and which converges very 
fast as well. A very interesting feature of this approach is 
its similarity with a density dependence of the pairing in- 
teraction. The numerical implementation of the present 
renormalization scheme is straightforward and amounts 
to very small changes of the existing codes. 
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